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Abstract 

Free Electron Laser (FEL) and Collective Atomic Recoil Laser (CARL) are described by the same 
model of classical equations for properly defined scaled variables. These equations are extended to 
the quantum domain describing the particle's motion by a Schrodinger equation coupled to a self- 
consistent radiation field. The model depends on a single collective parameter p which represents 
the maximum number of photons emitted per particle. We demonstrate that the classical model 
is recovered in the limit p ^ 1, in which the Wigner function associated to the Schrodinger 
equation obeys to the classical Vlasov equation. On the contrary, for p < 1, a new quantum 
regime is obtained in which both FELs and CARLs behave as a two-state system coupled to the 
self-consistent radiation field and described by Maxwell-Bloch equations. 
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I. INTRODUCTION 

Apparently very different systems as High-Gain Free Electron Laser (FEL) and Col- 
lective Atomic Recoil Laser (CARL) exhibit similar behaviors, showing self-bunching and 
exponential enhancement of the emitted radiation. Originally conceived in a semiclassical 
framework, they can be as well described quantum-mechanically . However, it 

is not explicitly evident how to obtain the classical limit starting from the quantum descrip- 
tion. First attempts to give a quantum description of the FEL have been proposed in the 
80s, starting from a canonical quantization of the A^-particle Hamiltonian in the Heisenberg 
picture, to study photon statistics and quantum initialization from vacuum in the linear 

n 



regime |^|9|. In 1988, Preparata proposed a quantum field theory of FEL [3], in which he 
has shown that, for ^ 1, the FEL dynamics is solved by a single-electron Schrodinger 
equation coupled to a self-consistent radiation mode. The same model has been recently 
obtained to describe CARL from a Bose-Einstein condensate (BEC) at zero temperature 
a a. Fu.he.n,o.e. it bee„ a.o p.oved expen.enta,ly Q Q Q that CARL . a 
BEC exhibits quantum recoil effects when the average recoil velocity remains less than the 
photon recoil limit. 

In this Letter we start from the classical model describing both CARLs and FELs and 
we extend it to the quantum realm showing the correspondence between the Preparata 
model and the CARL-BEC model 0. In particular, it is possible to derive an equation 
for the Wigner function of the A^-particle system. The Wigner function obeys to a finite 
difference equation which reduces to the classical Vlasov equation in the limit in which 
the number of photons emitted per particle is much larger than unity. In the opposite limit, 



both CARLs and FELs behave as a two-state system 

n 

Maxwell-Bloch equations |18l |. 
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described by the well-known 



II. CARL-FEL MODEL 



Apparently the physics of FEL and CARL appears to be quite different. The first de- 
scribes a relativistic high current electron beam with energy mc^ 70, injected in a magnet 
('wiggler') with a transverse, static magnetic field -B^ and periodicity A^, which radiates in 
the forward direction at the wavelength A ~ A^(l + a^)/27Q, where = eB^/mc^kw is 
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the wiggler parameter and = 27r/A^. Instead, CARL consists of a collection of two level 
atoms in a high-Q ring cavity driven by a far-detuned laser pump of frequency ujp which 
radiates at the frequency ~ cUp in the direction opposite to the pump. In both cases the 
radiation process arises from a collective instability which originates a symmetry breaking 
in the spatial distribution, i.e. a self-bunching of particles which group in regions smaller 
than the wavelength. 

It can be shown that, under suitable conditions and introducing proper dimensionless 
variables, the dynamics of both FELs and CARLs is described by the following Hamiltonian 

Q 

N 



n2 



- a^a, (1) 



where 6j and pj are the phase operator of the j-th particle and its conjugate momentum 
operator, obeying [Oj,Pj'] = iSjj'. In Eq.(PQ), a and are annihilation and creation operators 
for the forward radiation mode photon, with [a, a^] = 1. Notice that the dynamics described 
by Eq.(^ depends only on the parameter p and on the detuning 6, properly defined for the 
two systems: 

• For FELs, p = qpp and S = q'(7o— 7r)/7r, where q = mcfr/^k, 7^ = a/ (A^/2A)(1 + a^) 
is the resonant energy and pp = (l/7r)(aui/4cA;^)^/^(e^n/meo)^^'^ is the BPN parameter 

ri 

for a FEL 6 = [kyj + k)z — ckt, p = q{'-f — 7o)/7r and k = 2tt/X. 

• For CARLs, p = pc and 6 = {up — uj)/ujr, where ur = 2hk'^/m is the recoil frequency, 
PC = {So/uRY/^{u(Pn/2heoY'^, Sq = AQ/[2{r^ + + Q^)], Q is the pump Rabi 
frequency, A is the pump-atom detuning, F is natural decay constant of the atomic 
transition and d is the dipole matrix element Q]. Finally, 6 = 2kz andp = mVz/2hk, 
where Vz is the longitudinal atomic velocity. 

In these definitions, n = N/V is the particle density in the radiation volume V and m is the 
particle mass. Notice that in both cases p scales as n^/^, i.e. as the reciprocal of the inter- 
particle distance. Introducin^^j = {2/p)pj and A = {2/NpY^'^a, the Heisenberg equations 
associated with Eq.((H) are 



i=l ^ 

where j = 1 . . . and r = ujRpct for CARL whereas r = (47rpi?/A^)z for FEL. Considering 
the operators 6j, pj and A in Eqs.®-© as c- numbers, one obtain the well-known classical 
description for FEL and CARL. The scaling of Eqs.©-© is called 'universal' in the sense 
that, assuming resonance {i.e. 5 = 0), the equations do not contain any parameter. Hence, 
the scaling law of the various physical quantities can be obtained from their definition in 
term of p. In particular, from the definition of A and p it follows that the photon number 
per particle and the momentum recoil are proportional to p , both for FELs and CARLs. 
Hence, it is expected that when p ^ 1 the system behaves classically, whereas for p < 1 
quantum effects becomes relevant. Notice that \A\'^ + A^~^ 'YljPj ^ constant of motion in 
Eqs.(I21)-(|ll), i.e. the radiated intensity is due to the average recoil. 



III. QUANTUM CARL-FEL MODEL 



In ref. |3| Preparata, using quantum field theory, has shown that the collective dynamics 
of the system of A^ ^ 1 electrons in an FEL can be described by means of a single complex 
scalar quantum field whose behavior is governed by a Schrodinger-type equation in the 
self-consistent radiation field, which originates a pendulum-like potential: 

.dip 1 d'^ip W \ A id 



dr p Oe"^ 2 
dA f^^ i6 



^ [Ae^' - c.c] ip (5) 



, de\tlj\^e-'' + -A, (6) 
dr Jo p 

where ijj is normalized to one, i.e. J^^ d6\ijj{6,T)\'^ = 1. Note that Eq.Q is the Schrodinger 
equation associated to the Hamiltonian ([[]) and Eq.® corresponds to Eq.(@]) when the 
classical average of e"*^ is replaced by the quantum ensemble average. Quoting ref. 0|, 
Eqs.dni) and © are derived if one "formulate the many-electron problem in the language of 
quantum field theory and uses the large number A^ of electrons to evaluate the resulting path 
integral by saddle-point techniques". Rece ntly, the same model of Eqs.® and © has been 
used to describe CARL from a BEC llo| . Hence, we propose the nonlinear system of 

Eqs.dSI) and (jHl) as the quantum extension of the CARL-FEL classical model. We now show 
that the classical equations 0-0 are recovered in the limit p ^ 1. 
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IV. WIGNER FUNCTION APPROACH 



Let's consider the standard definition of tlie Wigner function for a state with wave func- 
tion tp{6, r): 

W{e,p,T) = ^ die'^^r [e - ^.r^i^ {o + \,r^ , (7) 

so that 

/ + 00 
dpW{e,p,T) = \^{e,T)\\ (8) 
-oo 

One can show that Eq.(j3I) is equivalent to the following finite difference equation for the 
quasi-probability distribution WiO.p.r): 

w {e,p+^,T^ -w {e,p-^,T^ =0.(9) 



dr de 2 

Using Eq. (jHl), Eq. (jH)) becomes: 

dA , ,« iS 



dr 



dp dew{e,p,T)e-'^ + -A. (10) 

■oo Jo P 

We underline again that Eqs.© and (fTUj) are equivalent to Eqs-dHl) and (jH)) using the Wigner 
function representation. In the right hand side of Eq. ©, the incremental ratio \W{6,p + 
e) - W{9, p-e)]/ (2e) dW{d, p) /dp when e = 1/p ^ 0. Hence, for p > 1 Eq. Q becomes 
the Vlasov equation: 

mV(e^p,r) SW(0p,r) _ gW(0 p,r) 

or o6 ^ ^ op 

Eqs. (fTUI) and (fTTj) are equivalent to the classical Eqs. ©-(El). This means that the particles 
behave classically, following a Newtonian motion, when p ^ 1 i.e. when the average number 
of photons scattered per particle is much larger than unity. In this limit, the quantum recoil 
effects due to the single photon scattering process is negligible. On the contrary, a quantum 
regime of CARL or PEL occurs when p < 1, in which each particle scatters only one photon. 
In fact, expanding the wave function in Fourier series as 

ij{e, r) = ^Y. Cn(r)e^'^(^+^\ n = -oo, . . . , +oo. (12) 

V 2,71 — 

and inserting this ansatz in Eqs.© and one can easily obtain the following closed set 
of equations for Qm,n{T) = Cm{r)*Cn{T): 

dQm,n _n)(S + m + n) g^.^ + ^ {.Qm+l,n - Qm,n~l) + A* (Pm,n+1 " Pm-l,n)|l3) 



OO 



dr p 2 
dA 
dr 




e 



FIG. 1: Numerical solution of Ea.()13|) and H14|) for p = 10 (first row), p = 1 (second row) and 
p = 0.2 (third row). The other parameters are 6 = 1, A{0) = 10~^ and c„(0) = 6no- Left column: 
dimensionless radiation intensity |Ap vs. r; central and right column: occupation probabilities 
Pn = |cnP vs. n and density distribution \'ip\'^ vs. 9, for r near the first maximum of \A\'^. 

where A = Ae"*^*^/^^"^. These equations are equivalent to Eqs.© and (jHI) for the density 
matrix in the momentum representation and have been discussed in ref. (l6|. FigQ shows 
the numerical solution of Eq.()13p and (fT^ for p = 10 (first row), p = 1 (second row) and 
p = 0.2 (third row). The other parameters are ^4(0) = 10~^, 0^(0) = 5„o and 6 = 1, which 
corresponds to a single photon scattering recoil. In the first column \A\'^ is plotted as a 
function of r. The central column shows the occupation probabilities Pn = |c„p vs. n, 
whereas the right column shows the density distribution lipl"^ vs. 6, for a value of r near 
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1(^A'S: + A'*S„) (16) 



the first maximum of We note tliat for p = 10 the system behaves classically: the 

momentum states are occupied in a range of the order of p and the average momentum is 
(p) ^ — p, as can be seen from the first row of figO Furthermore, the particle distribution 
shows periodic narrow peaks of density. When p = 1, the mainly occupied momentum states 
are those for n = and n = —1, corresponding to particles in the initial state or in the recoil 
state, respectively. Finally, when p <^ 1, the dynamics is that of a pure two-level system. In 
this limit, if the initially occupied state is the nth momentum state, the only two momentum 
states involved in the interaction are those for n and n — 1, so that Eq. p3p and (fT^ . after 
defining the 'polarization' Sn = 2^„_i „ and the 'population difference' Dn = Qn,n — Qn-i,n-i, 
reduce to the Maxwell-Bloch equations for a two-state system jisj ]: 

^ = -^A„5„ + A'D^ (15) 

dPn _ 

dr' 

where A„ = (5 — 1 + 2n)/p'^/^, A' = y/pA and r' = y/pr. With this new scaling and 
assuming resonance (i.e. A„ = 0), Eqs. (jl5|) - (fT7|) do not contain any parameter. Hence, the 
characteristic timescale is ruled by y/n instead of n^/^ as in the classical case. The quantum 
regime for CARLs and FELs is analog to the coherent spontaneous emission regime predicted 
quantum-mechanically for a two- level system in ref. 1^, where a series of optical "27r-pulses" 
are generated. In fact, assuming resonance (i.e. A„ = 0), A' and Sn are real. Hence, we 
can introduce the "Bloch angle" such that Sn = sin0, Dn = cos0. Then, Eqs. (|TH|) - (fT7jl 
reduce to a pendulum equation d'^(f)/dT''^ = sin(/) and d(f)/dT' = A'. Hence, in the quantum 
regime, the dynamics is that of a pendulum moving away from the unstable equilibrium 
point (0 = 0) and undergoing periodically a complete revolution ('27r-pulse') with angular 
velocity A'. 

Finally, we note that, adopting the same scaling of Eqs. (|TH|) - (|T7jl in Eq.ijH}, this can be 
interpreted as a Schrodinger equation for a single particle with a "mass" p'^^^ in a self- 
consistent pendulum potential. This provides an intuitive interpretation of the classical 
limit, that holds when the particle's 'mass' is large. The strong differences between the 
quantum and classical regimes are evident from Fig[3 
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V. CONCLUSIONS 



In this Letter we presented a unified quantum model that stands for apparently very 
different systems as FEL and CARL. The dynamics is described by a Schrodinger equation in 
a self-consistent pendulum potential and is ruled by an unique parameter p which represents 
the maximum number of photons scattered per particle and the maximum momentum recoil 
in units of the photon recoil momentum. The Schrodinger equation can be transformed in 
an exact equation for the Wigner quasi-probability distribution. The main results are the 
following: i) The classical model is recovered in the limit p ^ 1; this because the finite 
difference equation for the Wigner function reduces to the classical Vlasov equation, ii) 
In the limit p < 1 a completely different dynamical regime occurs (see Fig. Q): due to 
momentum quantization the system reduces to only two momentum states obeying to the 
Maxwell-Bloch equations which describe the dynamics of a two-level atomic system coupled 
to a coherent field. 
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